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PREFACE 

This final report documents a study which examined the scientific 
evidence for the addition of a 7th spectral band to the Landsat-D 
Thematic Mapper (TM) . This study was performed within the Infrared and 
Optics Laboratory of the Environmental Research Institute of Michigan. 

The research covered in this report was performed under Task Order No. 

One of NASA Contract NAS9-15362 during the period from August - December, 
1977. Mr. Robert K. Stewart o£ the Earth Resources Program Office at 
Lyndon B. Johnson Space Center, Houston, Texas was technical monitor 
for this task. 

The work performed under this task represents a "quick response" 
effort in which several candidate bands proposed by members of a large 
"user community" were examined from the standpoint of the widespread 
usefulness of the resulting imagery from a spaceborne sensor . This 
study would have fallen far short of the intended goal without the 
cooperation and sharing of technical data from several research groups 
made possible by the good offices of the technical monitor and the Land- 
sat-D Project Scientist Dr. Vincent Salomonson. Even then, the lack 
of adequate spectral measurements in certain spectral regions prevented 
a comparable chance from being given to each of the candidate spectral 
bands . 

This task was conducted under the guidance of Mr. Richard R. Legault 
Head, Infrared & Optics Division. Dr. Quentin Holmes was the project 
manager for this study. The authors wish to acknowledge the substantial 
contributions made to this work by Messrs.: Jerry Beard, Richard Nalepka, 
Robert Horvath, Bill Smith, Lester Witter, and George Lindquist of ERIM. 

We wish also to thank Dr. Alex Goetz and Dr. Anne Kahle of NASA/JPL for 
making available to us the excellent field measurement data they have 
laboriously gathered and analyzed over a period of several years, and 
Dr. Robert Vincent of Geospectra Corporation for his consultation in the 
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use of the thermal infrared region. Dr. Stephen Ungar and Ms. Vivien 
Gornitz of the Goddard Institute for Space Studies shared preliminary 
results from their studies with us, while Dr. John Barker and Dr. Mel 
Podwysocki briefed us on related activities in their groups. Thanks 
also go to Dr. Larry Rowan/USGS and Dr. Fred Henderson from the GEOSAT 
committee. Finally, thanks must also be given, to the many researchers 
who took time out from their busy schedules to respond professionally 
to our telephone survey in order that the significant issues involved 
in the selection on a seventh band could be identified. 
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1 

EXECUTIVE SUMMARY 

In the planning for the Landsat-D Thematic Mapper provision was 
made for the incorporation of a 7th spectral band provided the increase 
in technical performance would justify the increased cost- The variety 
of possible applications combined with the diversity of professional 
opinions within a given application and the intense interest of a large 
community of users of spaceborne imagery required that an assessment 
be made of the underlying scientific evidence before a final recommen- 
dation was put forth. 

Beginning with a review of existing ideas and an extensive litera- 
ture survey by a senior researcher, a telephone survey was made of more 
than thirty well-known individuals who are currently active in remote 
sensing research to ascertain the critical issues in their discipline 
which needed to be considered in selecting a 7th band for the Thematic 
Mapper. The diversity of professional belief even within a specific 
application appears to be a basic characteristic of remote sensing 
technology. In order to ensure maximum information flow, the ground 

A 

rule was adopted, that a respondent's name would not be released with 
their comments. To permit the validity of conflicting concepts to be 
assessed each respondent was asked to provide primary sources (based 
upon their current work whenever possible) which would substantiate 

key points. Based upon a consensus of the people interviewed, the follow- 
ing proposed spectral bands were examined in detail: 

- 0.70-0.7ii micrometer spectral band for agricultural applications 

- 1.20 micrometer centered spectral band for forestry & agricul- 
tural applications 

- 2.20 micrometer centered spectral band for geological applications 

- 8. 2-9. 2 micrometer spectral band ("2nd thermal band") for geolog- 
ical applications. 
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Each of these were examxned in terms of the feasibility of gath- 
ering high quality imagery from space while taking into account solar 
illumination, atmospheric attenuation, and the expected signal/noise 
ratio achievable within the sensor constraints of the Thematic Mapper 

design. For the 2.2 ym band and the 2nd thermal band the availability of 
field and laboratory spectral measurements permitted the performance 
gains to be quantified in terms of the increase in percent correct clas- 
sifications of hydro thermally altered rocks and broad classes of rock 
types, respectively. The technique employed to compute percent correct 
classification was a linear discriminant analysis applied to inband 
signal values. The inband signal values were calculated from represent- 
ative spectral reflectance/emittance curves for specific samples. Solar 
illumination, atmospheric absorption and attenuation as well as sensor 
noise characteristics were taken into account by the radiative transfer 
model used to transform ground based reflectance values into the cor- 
responding sensor signals recorded onboard a spacecraft. 

We emphasize that, unfortunately, no comparable test could be 
given to all the candidate bands which were suggested . The absence of 
sufficient reflectance/emittance spectra forced us to this proceeding. 
Based upon analysis of the state of the art at the end of 1977 the 
following conclusions are drawn and recommendations are made: 

1. Although technically feasible, a band in the 0.70-0.74 micro- 
meter region is not recommended for incorporation on the Land- 
sat-D Thematic Mapper because, at present, there is insuffi- 
cient evidence to justify its inclusion. 

2. Evidence for incorporating a band in the 1.2 micrometer region 
(there are water absorption bands at 1.15 and 1.45 micrometers) 
is inferential and rests upon the idea that having two 
observations in the near infrared plateau would permit a 
resolution element to be better characterized in terms of a 

2 " 
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mixture of vegetation and soil on different soil types. This 
band was no£ selecte'd because low detector sensitivity com- 
bined with the already high performance of the TM in agricul- 
tural applications made this spectral band a 3rd choice. 

A single thermal b'and on the Landsat-D Thematic Mapper will 
probably not provide the large improvement in results which 
is commonly expected because of the competing physical phenom- 
ena. 

Incorporation of a 2nd ’'thermal band (8.2— 9.2 micrometers) 
would improve classification accuracies for some general rock 
types. However, using the limited amount of emission spectra 
available, our work does not support the contention that accur- 
ate identification of rocks according to their S 1 O 2 content 


would be possible. 

- In order to simplify the data processing (make the data more 
accessible to the user) , a 90m ground resolution for the ther- 
mal band is preferable. Unfortunately, the selected numbers 
of detectors cannot be changed at this late date. Great care 
will be required in the way these thermal data are used in 
conjunction with the reflective bands during machine process- 
ing. 

- It is recommended that a program be undertaken to collect 
field and laboratory spectral data spanning both the reflec- 
tive and emissive portion of the spectrum so that in future 
studies such as this the number and location of spectral bands 
in the thermal region needed to do accurate identification of 
general rock types can be quantified. 

- Because the performance of a 2nd thermal band for general geo- 
logic mapping is on the same order as that obtainable with a 
2.2 ym-^band, this candidate spectral band is considered to be 
a 2nd choice. 
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4. Based upon the substantial improvement in mapping zones of 

hydrothermally altered rocks from unaltered zones under a broad 
range of illumination and atmospheric conditions (from 78% 
accuracy with the original TM bands to 92% accuracy with thxs 
band) a 2.08-2,35 micrometer spectral band havxng .a ground 
resolution pf 30 meters is recommended for inclusion on the 
Landsat-D Thematic Mapper. Imagery from this spectral band will 
also be helpful for general geological mapping, as well as for 
discrimination among various classes of vegetation with different 
degrees of vigor. 


4 
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2 

BACKGROUND AND INTRODUCTION 

The current set of six Thematic Mapper spectral bands are the 

JL 

end result of studies and discussions which have taken place over a 
period of several years. One of the most significant of these activi- 
ties was the Landsat-D Thematic Mapper technical working group convened 
at NASA/Johnson Space Center in mid-1975. Four panels, each composed 
of a number of remote sensing specialists with a variety of backgrounds 
were asked to independently develop recommended specifications for the 
Thematic Mapper with the understanding that the primary application 
for this sensor would be agriculture. Three out of four of these 
panels came in with recommendations for the allocation of two spec- 
tral bands to the near infrared reflection plateau (0.74-1.3 micro- 
meters) of vegetation. The remaining spectral bands were allocated 
to the green and red portion of the visible region, the atmospheric 
window near 1.65 micrometers and the thermal infrared region around 
11 micrometers. Based upon their understanding of the needs of poten- 
tial users a spatial resolution of 30-40 meters ground IFOV combined 
with a radiometric accuracy on the order of 0.5% NEAp (noise equiva- 
lent change in ground reflectance) was recommended for all of the 
reflective bands except the 1.65 micrometer band where only a 1% 

NEAp was achievable. For the thermal band the recommended radiometric 
sensitivity was 0.5°K NEAT (noise equivalent change in ground tempera- 
ture) although It was clear that this would probably only be achieved 
at the expense of accepting a larger ground IFOV for this band than 
for the other spectral bands. 

In subsequent work the precise locations and spectral bandwidths 
were chosen in such a manner than the two bands allocated to the near 
Infrared region were placed side by side (0.74-0.80 and 0.80-0.91 micro- 
meters, respectively). Not surprisingly, these two bands were 

*EOSPDG, 1973, Thompson, F. et al; 1974, Hone3!well Radiation Center; 

1974, Harnage, J.; 1975, Mallla, W. et al; 1976, Morgenstern, J. ; 

1976, CORSPERS. 
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shown to have extremely high cp.rrelation (-R = 0 . 99 ) over all the data 

sets developed in an ERIM study (Malila, et al, 1976) which sought to 
quantify the expected performance of the Thematic Mapper through sxmula- 
txons based upon high quality aircraft imagery. Because two spectral 
bands wxth such extreme correlatxons are redundant, they were later com- 
bxned and a new spectral band xn the blue region was chosen. At that 
time many remote sensing specialxsts regarded thxs new spectral band 
(0.45-0.52 micrometers) as useful prxmarxly for hydrological applications 
and research. However, a detaxled study (Morgenstern et al, 1976) xn the 
Spring of 1976 showed that to achieve high classif xcation accuracies with 
maximum likelxhood processing over agricultural scenes each of the five 
reflectxve bands on the Thematic Mapper would be important at some txrae 
durxng the growing season. In that study, TM band 1 was found to be more 
xmportant than either TM band 3 or the 120 meter resolutxon thermal band 
when identifying corn and soybeans in mid- July. 

In the Fall of 1976 the National Research Council's Committee on 
Remote Sensing Programs for Earth Resources Surveys (CORSPERS) completed 
an xndepth review of the Thematic Mapper as the follow-on sensor to the 
present multispectral scanner carried by the Landsat spacecraft. Based 
upon thexr assessment of the results of remote sensxng research and the 
underlyxng science xn applxcatxons ranging from Cartography, Inventorying 
and Assessment of Vegetatxon, and Land Use Management to Geology, Ocea- 
nography, Water Resources Management, and Environmental Monitoring, thxs 
group's recommendations (CORSPERS, 1976) were for only slight changes in 
the location and band widths of the reflective bands of the Thematxc 
Mapper. However, in the thermal region they recommended a broader ther- 
mal band (8.8-12.6 ym as opposed to NASA's choice of 10.4-12.5 ym) to per- 
mit use of a 90 m ground IFOV for this band in place of the 120 meters 
selected by NASA. 

Figure 1 lists the spectral bands for the Thematic Mapper whxch 
were put forth by NASA shortly after the Landsat-D technical working 
group xn 1975, the bands recommended by CORSPERS in the Fall of 1976, 
and the spectral bands currently specif xed for this sensor. 
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Landsat-D Thematic Mapper 
Technical Working Group 
(mid-1975) (urn) 

Early 1976 
(urn) 

CORSPERS 
(Fall 1976) 
(ym) 

Current Thematic 
Mapper Sensor 
Specification* 
(ym) 


0.45-0.52 

0.47-0.52 

0.45-0.52 

0.52-0.60 

0.52-0.60 

0.53-0.58 

0.52-0.60 

0.63-0.69 

0.63-0.69 

0.62-0.68 

0.63-0.69 

0.74-0.80 ^ 

0.80-0.91 ) 

0.76-0.90 

0.76-0.90 

0.76-0.90 

1.55-1.75 

1.55-1.75 

1.55-1.75 

1.55-1.75 

10.40-12.50 

10.40-12.50 

8.80-12.60** 

10.40-12.50 

ic 

As of August 1977- 









Chosen to reduce ground IFOV of the thermal band from 120 meters to 
90 meters . 


FIGURE 1. SPECTRAL BANDS RECO MNEND SD FOR THE L&EDSAT-D 
THEMATIC MAPPER BY VARIOUS GROUPS 


For completeness 5 we now summarize the major conclusions prior to 
the current study: 

Initially, the Thematic Mapper spectral bands were selected with 
the agricultural application in mind. 

- A spatial resolution of 30 meters ground IFOV was selected in 
order to overcome some of the limitations of the current Landsat 
MSS without incurring the difficult sensor design and large data 
processing loads associated with an even smaller IFOV. 

In an ERIM study radiometric sensitivities poorer than; 

NEAp ~0.5% for the first four spectral bands 
NEAp ~1.0% for the 1.65 micrometer spectral band 
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and NEAT ~0-5°K for the thermal- band were shown to result in 

a significant degradation of maximum likelihood classifica- 
tion accuracy over agricultural regions . 

- With the exception of the poor spatial resolution of the thermal 
band, CORSPERS found the Thematic Mapper spectral bands closely 
paralleled the bands they would recommend for a variety of appli- 
cations. However, they too had a long list of recommendations 
for additional spectral bands. 

The purpose of the current study is to survey the evidence which 
supports several candidates for a 7th band and quantify the performance 
gains which would be achieved if a 7th spectral band were added to the 
Landsat-D Thematic Mapper. 


8 
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3 

APPROACH 

Since the Landsat-D Thematic Mapper may well fulfill a role in 
-the 1980 's quite similar to that played by the current Landsat MSS in 
the 1970 's it is important that recommendations for any additional 
spectral band be based upon the best available scientific evidence. 
Although we began by summarizing the existing literature, the need 
was to go beyond the level of information available in the published 
literature in order to assess the significance of a potential appli- 
cation, and to identify the primary sources of technical data which 
supported the selection of a spectral band for this application. To 
us this required that a frank discussion be held with a number of 
individual scientists who are actively engaged in remote sensing 
research. The intent was to ascertain their professional opinions 
and concerns regarding the critical issues involved in the selection 
of a 7th band for the Thematic Mapper. 

For each candidate band proposed the underlying concepts and 
evidence were examined, and from the large number of candidates a sub- 
set was identified for even more detailed analysis. 

For this subset of candidate bands calculations were made of the 
noise equivalent change in radiance which could be expected under a 
variety of sun angle and atmospheric conditions taking into account 
state-of-the-art detectors and assuming current Thematic Mapper values 
for electronic bandwidth, number of detectors /channel, and orbital 
parameters. The results of these computations were then compared to 
the magnitude of the change in radiance, at spacecraft altitude, pro- 
duced by the phenomena being detected. 
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For candidate bands in two of the spectral regions, representative 
laboratory and field spectra were used in conjunction with a radiative 
transfer model to produce inband radiance values under the wide range 
of conditions of observation for which this sensor will be utilized. 
These values were then digitized taking sensor noise into account and 
a linear discriminant analysis was used to quantify the improvement 
in classification accuracy obtained by addition of a candidate spec- 
tral band. 


10 
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4 

SURVEY OF REPRESENTATIVE RSIOTE SENSING SPECIALISTS 

As presented in the introduction, a, senior ERIM researcher 
contacted a representative sample of scientists who are actively 
involved in remote sensing research to ascertain their professional 
opinions and- concerns regarding the applications to be pursued, the 
appropriate spectral regions, and the existence of technical data which 
would support their suggestions. In order to permit open information 
flow, we adopted the ground rule that no responses would be released 
with the contributor's name attached. Some of the opinions obtained 
were in direct conflict with official agency policies and had to be 
handled as privileged information. This is an unusual proceeding and 
illustrates the tension associated with the scientific and political 
importance of a 7th band Thematic Mapper decision. 

The different professional opinions about the usefulness of 
adding a specific candidate band to the Landsat-D Thematic Mapper 
manifested itself in comments ranging^from high consent to total 
rejection. The different interests have clashed repeatedly to the 
extent that NASA's plan to extend the flow of orbital earth observa- 
tion data by using the Thematic Mapper to upgrade the quality of that 
data has not been accepted by all scientists. Concerned about relying 
totally on this new sensor, they wish to maintain an MSS on Landsat-D 
and support a design to meet operational requirements. 

With the selection of inventorying and assessment of vegetation 
as the primary optimization objective in mission design, the program 
ran the risk of not serving equally the large Landsat user community 
which has developed. So, after having taken into consideration a 
band between 0.45 and 0.52 micrometers (ym) for water research, it 
seems reasonable that the addition of a spectral band for geological 
applications would be suggested. 
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Even in the field of geology, two main opinions crystallized: 

^ - One Supporting a band around 2.2 ym, based upon the 

research of L. Rowan and A. Goetz, etc. 

- The other supporting a second Thermal IR channel in the 
8. 0-9.0 ym wavelength region, based upon the research of 
R. Vincent and R. Lyon. 

Figure 2 is a synthesis of professional opinions by cognizant 
remote sensing scientists about the possible selections of a 2.0-2. 6 ym 
spectral band for geological remote sensing and what might be preferred 
as an alternative. 

The consensus of these contacts with the different scientists, 
was that the 2. 0-2. 6 ym wavelength region provides information on 
altered zones. This application is not significant to all geologists, 
but very important to economic geologist looking for metals. Some 
of the geologists, who are more interested in lithology and structure 
feel hydrothermal alteration is a limited application, not specific to 
the Office of International Geology or developing country needs, and 
would prefer the thermal IR or even the current MSS. It was also pointed 
out quite clearly, that the worst eventuality for the geology community 
would be, that neither the 2. 0-2. 6 ym band nor the second thermal IR 
band will be added to the Thematic Mapper. Instead of debating about 
the better band the ideal solution would be to fly both, a 2.0-2. 6 ym 
and a second thermal IR band (8. 2-9-2 ym) . This of course could only 
be accommodated in the current design by eliminating an existing band 
which would simply cause other problems. 
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PRO 


CON 


2.2iain centered band 


-Enthusiastic about 2.0-2. 6 litn 
because of Rowan's and Goetz's 
work. 

-Very important for separating 
chemically altered rocks 
-JPL's work is convincing 
-Won't produce apples and or- 
anges . 


-Not enthusiastic about 2. 0-2. 6 ym, 
but better than nothing 

-Very limited for geological mapping 
use. 

-Results are not convincing 

-Too broad, will produce apples and 
oranges . 


00 


Second Thermal IR band -Prefers two thermal bands 


-Very important for quartz/ 
calcite discrimination, silica 
content, provides redundancy for 
thermal information and tempera- 
ture measurements. 


-Decision on Thermal IR should wait 
until after we can study Landsat-C 
experiment . 

-Not enough experience, more research 
needed before optimum IR band can 
be selected. 


NOTE: There were also a few scientists who believed that other wavelength regions and 

applications were more important such as : 

- Replacement of the current 1.55-1.75 pm band by a 1.55-1.65 ym and a 1.65-1.75 ym 
band. This is most important to identify geological alterations. 

- 2.70-3.0 ym band which provides information on quartz absorption and better dis- 
crimination of sedimentary rocks. 

- Most interested in 3. 0-3. 5 ym or 3.0-3.15 ym region for daytime and nighttime 
comparisons. 

- Prefer photography to either of the above. Both spectral regions need further 
research for geological applications. 


FIGURE 2. SYNTHESIS OF PROFESSIONAL OPINIONS CONCERNING 
A 7TH SPECTRAL BAND ON THE LANDSAT-D THEMATIC MPPER 
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5 

CLOSING THE GAP BETTOEN PUBLISHED DATA Ai^D CURRENT KNOWLEDGE 


To' base th'e recommendation for a 7th band on the best technical 
information availaljie at this time, the existing literature was supple- 
mented through meetings with' groups which are currently conducting re- 
lated research at NASA/JPL, and NASA/GSEC. Arrangements were also made 
for use of as yet unpublished technical data from the following sources: 
Earth Resources Spectral Information System (ERSIS) at ERIM 

- Chevron Oil Corporation 

- NASA/ Jet Propulsion Laboratory (JPL) 

- Geospectra Corporation 

— Goddard Institute of Space Studies (GISS). 

Key concepts contributed by each of these groups will be presented in 
subsequent sections where each candidate is treated. Specific data 
obtained from these sources was as follows: 

- NASA/JPL provided us with a large, representative, sample of 
calibrated rock spectra in the 0.4-2. 6 pm region. These spectra 
provided an excellent experimental basis for the quantitative 
analysis of the performance of a 2.2 pm centered band. 

- Chevron Oil Corporation released rock spectra and analysis results 
obtained from processing M-7 aircraft imagery from the Cady Hills 
region of California. These results were instrumental in bring- 
ing about the re-examination (see Appendix A) of the underlying 
physical processes responsible for spectral structure in the 

2-3 pm region. 

- Geospectra contributed several, previously unpublished, thermal 
spectra which were used in support of the thermal studies. 

- Goddard Institute of Space Studies made available the agricultural 
field measurement data in the 0,4— 1.0 pm region used to support 
the 0.70-0.74 pm candidate band. They also shared with us the pre- 
liminary results obtained from digital processing of 24 channel 
scanner data collected over Waterpocket Fold, Utah. 
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Goddard. Space Flight Center briefed us on' their results using 
24 channel scanner data over a number of geological sites and on 
the spectral band selection technique they have been developing. 


16 
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RADIATIVE TRAHSEER CALCULATIONS 

In order to compute the changes in radiant energy at. spacecraft 
altitude which were associated with the discrimination of specific tar- 
gets It was necessary to employ a radiative transfer model which took 
into account aerosol scattering and molecular absorption as well as solar 
alltmination and viewing, geometry . The model employed in making all of 
the computations in the current study was an extension of Turner's Radi- 
ative Transfer Model such that it provided a consistent model which is 
applicable to a wavelength region from 0.4 to 15 micrometers. The 
original model (Turner, 1974) was modified to include ozone absorption 
in the vicinity of 0.5 micrometers. Then, it was combined with the Aggre- 
gate Model for molecular absorption.- In the process it was necessary to 
make some simplifying assumptions in order to permit the numerical evalu- 
ation of the integrals involved. 

Heuristically, we note that due to strong forward scattering (as 
evidenced by the, dominant peaking of the scattering phase functions in 
the forward direction) : . 



the scattered path length the direct path length. 


Hence, the solar irradiance extinction can be represented as a product 
of the attenuations due to atmospheric scattering and molecular absorp- 
tion. 


T = T 

total scattering 


absorption 


( 1 ) 


Moreover, under the additional assumption that the altitude distribution 
of scatterers and absorbers is identical, path radiance values can be com- 
puted directly. This last assumption is quite reasonable for haze levels 
corresponding to visual ranges greater than 13 km. This assumption is 
not so well taken for severe haze levels (i.e. , visual ranges < 5 km) 
where the majority of the scatterers are concentrated in the lower por- 
tion of the atmosphere. 
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Figure 3 shows a typical curve of spectral radxance in units of 
-1 -2 -1 -1 

photons sec cm sr Um at the top of the atmosphere which was 
calculated 'using this model. This particular curve is based upon the 
reflection/ emission sp'ectrum 'of a healthy deciduous tree whose tem- 
perature is 300“K. The solar zenith angle for these computations was 
45° and the sensor zenith angle is 0° (i.e., looking at nadir) the atmos- 
pheric composition is characteristic of a tropical region during midday 
with a haze level corresponding to a visual range of 23 km (i.e., a 
reasonably "clear day") . A logarithmic scale was chosen for the vertical 
coordinate 'of this plot ’because irradiance from the sun peaks around 
0,6 urn and falls off very rapidly toward longer wavelengths. 

The "green peak" at 0.55 ym, the chlorophyll absorption dip in the 
red region -from 0.67 to 0.68 ym and the high near IR plateau are all 
characteristic of a vegetative target. While the absorption dips at 
1.15 ym, 1.45 ym, and 1.95 ym are due to atmospheric water, the dip at 
1.38 ym is due to CO^. The atmosphere is essentially opaque in the 
region from 2. 5-2. 8 ym due to a number of overl apping H^O and CO^ absorp- 
tion bands. In the emissive region between 3 ym and 14 ym the envelope 
of the spectral radiance is essentially that due to emission by a gray 
body with an emissivity of 0.95 at 300 °K with a number of atmospheric 
absorption bands superimposed on it. At 300°K spectral radiance due to 
emission from the target rises to a broad maximum at around 10 ym. In 
addition to apparent temperature differences some materials have reststrah- 
len bands which serve as a basis for discrimination in the thermal region. 

In the reflective region of the spectrum our computations of inband 
radiance were carried out by first calculating the spectral radiance at 
the top of the atmosphere for two limiting situations, 100% ground reflec- 
tance and 0% ground reflectance. The total radiance, L(A) , under a 
specified condition of observation when viewing a material of arbitrary 
spectral reflectance, p(A), is then given by. 



RflOIRNCE (PH0 TQNS/SEC/CMhk2/3TER/UM) 

rxiois Mtoi‘1 ixiois ixi&te _ ixio” 


TRQPICf!L(HIDDftY) VIS RRUGE «• 23 KM, VIEW ZEM flNG - 


0.0 QEG, SOLflR ZEN flNG « 


tiS.O DEG, TOr TEMP - 


300. 0 K. n 



FIGURE 3. A TYPICAL CURVE OF THE RAMAUCE AT THE TOP OF THE ATMOSPHERE PREDICTED 
BY THE RADIATIVE TRANSFER MODEL. Thxs axample was computed from 
the spectra of a healthy decxduous tree. 
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L(X) - (X) - (X)] , p(X) + t^(X)} (2) 

where (A) and (A) denote the total spectral radiance at the top 

of the atmosphere for a 100% and a 0% reflecting surface under the speci- 
fic conditions of observation (i.e., represents the path radiance). 

The factor hc/A (where h and c are Planck's constant and the velocity of 
light, respectively) converts the units from photons /sec into watts. Thus, 
the spectral reflectance, p(A), for a given material is transformed into 
the corresponding total radiance at the top of the atmosphere by a two 
parameter family of functions. This is illustrated in Figure 4. 


In practice, both t^^^CA) and t^(^) are computed by numerical integration 

once and for all using the radiative transfer model for a specific at- 
mospheric constitution, haze level, sun zenith angle, and sensor viewing 
geometry. Then Eq. (2) can be used to efficiently compute spectral 
radiance, L(A), under these same conditions for any material of whose 
spectral reflectance, p(A), is known. 

The counterpart of Eq. (2) in the emissive region is given by, 

L(;,) = ^ ’^path^^^} 

where L(t,A) is the radiance emitted by the material being viewed on the 
earth's surface. This equation, which reduces to Eq. (1) when L(t,A) is 
negligible, is a linear equation in three parameters: solar radiance 

[L^(A)J, path radiance 1 radiance from the ground [L(x,A)]. 

Hence, computations of inband spectral radiance in the thermal region 
must be made with the radiative transfer model for each individual 
material . 
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FIGURE 4. RADIATIVE TRANSFER MODEL GEOMETRY 
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ESTIMATION OE SENSOR NOISE LEVELS 


One of the considerations involved in assessing the suitability of 
a given band for a spaceborne scanner is its noise equivalent signal 
level. In the course of this study estimates of the noise levels which 
could be expected were made for a variety of spectral bands . In all 
instances we assumed that the detectors employed would be individually 
selected, state-of-the-art devices. To allow for uncertainties in 
estimating the achievable values for key parameters such as transmission 
of the sensor optics, and detectivity of the detector, computations 
were made for both a pessimistic case and an optimistic case. 

A straightforward development of the physics involved ( See Appen- 
dix leads us to the well-known expression for noise equivalent power, 
NEP , as 

NEP . (4) 


where A is the area of the detector, and B is the electrical band- 
d 

width of the sensor and D* is the detectivity of the detector. Using 
this equation for noise equivalent power the corresponding change in 
ground reflectance can be expressed in terms of the atmospheric transmis- 
sion [t^], the sample irradiance [H ] , transmission of the optics [t ], 
and the f-number resulting expression is 


NEAp = 


4F 


T.H^T D' 
A I O 



(5) 


Equations (4) and (5) were used in conjunction with the Landsat-D 
Thematic Mapper parameters shown in Figure 5 to estimate sensor noise 
levels. Representative values for D* were chosen based upon an assess- 
ment of the state-of-the-art of detector technology by the Santa Barbara 
Research Center (Bode, 1976). 
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Orbital Altitude 
Orbital Velocity 
Swath Width 
IFOV 

Detectors/ channel 


705 km 

7,615 meters/sec 
185 km 

30 meters (120 meters for the thermal band) 
Channels 1-5:16 detectors, channel 6: 4 detectors 


Min Radiance (S/N) Max Radiance (S/N) 


Spectral 

band 

1 

0.28 mW/cm^/sr 

(32) 

1.00 mW/cm^/sr 

(85) 

Spectral 

band 

2 

0.24 mW/cm^/sr 

(35) 

2.33 mW/cm^/sr 

(170) 

Spectral 

band 

3 

0.13 mW/cm^/sr 

(26) 

1.35 mW/cm^/sr 

(143) 

Spectral 

band 

4 

0.16 mW/cm^/sr 

(32) 

3.00 mW/cm^/sr 

(240) 

Spectral 

band 

5 

0.08 mW/cm^/sr 

(13) 

0.60 mW/cm^/sr 

(75) 

Spectral 

band 

6 

300“K 


320“K 



FIGURE 5. NOMINAL PARAMETERS OF THE LANDSAT-D THEMATIC MAPPER 


The technical difficulties associated with the machine processing 
of scanner imagery which has a 120 meter IFOV for the thermal band and 
30 meter IFOV for the bands in the reflective portion lead the authors 
of this report to inquire if the number of detectors used for the 
thermal band could be increased to provide say 90 meter IFOV. Unfor- 
tunately, the selected number of detectors for the reflective region 
is not divisible by three, and the time for the major revision in the 
sensor design required to make this possible has long since passed. 
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CANDIDATE BAWD IN THE 0.70-0.74 ym REGION 

8 . 1 Background 

Based upon a complex, mathematxcal investigatxon of scene charac- 
terxstics by a group of Soviet scxentists (Kondratyev, K. et al, 1975) 
in which the 0.70-0.74 ym region was xndicated as being hxghly signifi- 
cant, a band in thxs regxon was consxdered as a candxdate 7th band for 
the Thematic Mapper. A related, but considerably easxer to interpret, 
development by an*Amerxcan group of scientxsts at Goddard Instxtute of 
Space Studxes appears to support this same regxon for agricultural appli- 
cations involving the discrimination of wheat from other agrxcultural 
crops . 

8.2 Avaxlable Evidence 

Seeking to xdentify the spectral regions most effective in charac- 
terizxng remotely sensed scenes, Sovxet researchers applxed the mathe- 
matxcs of informatxon theory to the results of remotely sensed experi- 
ments from a broad spectrum of different applicatxons . In vxew of the 
common practxce of sensor designers to avoxd the shoulders of reflection 
peaks, many researchers were surprised when the results of thxs study 
xndicated the 0.70-0.74 ym region should be intentionally included in a 
sensor designed to be used for agricultural applications. Even after a 
careful review of this work as it is reported in the literature, we can 
not offer a simple physical explanation as to why the 0.70-0.74 ym band 
was emphasized. 

Based upon an interpretation of high resolution spectral radiance 
measurements collected over the central portions of agricultural fields 
in the Imperial Valley of California (Ungar, S. et al, 1977), researchers 
at the Goddard Institute for Space Studies believe that a narrow spectral 
band i,n this region would permit wheat in the heading stage to be distin- 
guished from a number of other agricultural crops. Typical curves used 
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to support thxs concept in the briefing given to us by Dr. Stephen G. 
Ungar are shoxm in Figure 6. The curve labeled "Wheat" was obtained 
by averaging together spectral radiance curves measured while observing 
the interior of a single wheat field in the headed stage with a modified 
Jarrell- Ash 0.33 meter spectrograph mounted in a light aircraft. These 
"field measurements" provide digital value of radiance in 410 channels 
each having a bandwidth of 0.0014 pm in the spectral region from 0.43-1.0 
micrometers. The difference between the spectral radiance obtained 
when viewing headed wheat and when viewing alfalfa is indicated by the 
arrows labeled "far red shift". Based upon his research, Ungar recom- 
mends using a band as narrow as possible, and no broader than 0.2 ym. 

For the curves shown, the largest differences occur over the region from 

0.71-0.73 ym. Numerically integrating these differences over this 

-2 -1 

spectral range gives a value of about 40 microwatts cm sr . Although 
this change in radiance levels is more than sufficient to be measurable 
at spacecraft altitudes, the results are inconclusive. The fact that 
these measurements of observed spectral radiance were never converted 
into corresponding values of spectral reflectance precludes comparisons 
with other curves in this atlas which were collected under different 
conditions of illumination. Thus, we were unable to assess the consis- 
tency of this interpretation for even the crop types which were observed 
in the Imperial Valley. 

While there are differences of professional opinion (Private commu- 
nication from C. J. Tucker, October 7, 1977) about the existence and 
applicability of this "far red shift" phenomena, we note that a 0.71- 
0.73 ym spectral band was flown on the ERIM M7 airborne scanner back in 
1972, Based upon notions similar to those investigated by Stephen Ungar, 
this spectral band was included as a "wild card" on a flight over an 
agricultural region in Canada. Unfortunately, time and resources did 
not permit even the limited analysis to be carried out which was 
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originally envisioned by the researcher who suggested this band. 

Although collected from an aircraft, this scanner imagery has signal/ 
noise values comparable to data which would be collected from satellite 
altitude with the Landsat-D Thematic Mapper. Since this imagery existed, 
it was examined in the course of the present study. Figure 7 shows two 
spectral bands (0.71-0.73 ym and 0.55-0.60 ym) collected by the M7 
scanner over Guelph, Ontario" in July "1972 . As expected^ the 0.71- 
0.-73 ym band does not exhibit as much contrast as a spectral band which 
corresponds closely to TM band 2. When the appearances of several agri- 
cultural fields in the 0.71-0.73 ym spectral band were visually compared 
with their appearance in the 1.55-1.80 ym band and the 2. 0-2. 6 ym band 
it was clear that this band contained different information. Unfortun- 
ately, the ground truth which was collected for this flight was not 
suitable for use in quantifying this additional information content. 

8.3 Conclusion 

Although technically feasible, a spectral band in the 0.70-0.74 ym 
region is not recommended for incorporation on the Landsat-D Thematic 
Mapper because, at present, there is insufficient evidence to justify 
its inclusion. 
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(a) 0.71-0.73 ym 



(b) 0.55-0.60 ym 


FIGURE 7. ELLORA RESEARCH FARM, GUELPH, ONTARIO, CANADA 
ACQUISITION DATE: 6-17-1972, ERIM M7 SCANNER 
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variabilxty of spectral radiance values in this region due to different 
concentrations of atmospheric water can amount to as much as 30% of the 
nominal signals). Dr. Hoffer then pointed out that the result from this 
study alone was certainly not adequate to justify a conclusion. However, 
It did serve as another important piece in a puzzle. In his opinion, 
this result, in conjunction with similar results obtained during other 
experiments, seemed to indicate the importance of having more than one 
spectral band in the near IR region where vegetation has a relatively 
constant reflectance. More recent work at LARS (Bauer, M. E. et al, 1977) 
has suggested that an additional spectral band in the 1.20-1.30vim region 
would be useful in mapping soils in the presence of vegetation. 

PbS detectors provide the highest values of detectivity in the spec- 
tral region of 1.2pm. For selected detectors it would be possible to 

10 

obtain D* values of 6 x 10 in an uncooled environment. Unfortunately, 
PbS detectors are not usable for this particular application because 
their response time is too slow. To provide a 1.2-1.3hm spectral band 
with a 30 meter IFOV using 16 detectors/band the need is for a detector 
with a bandwidth of 51.2 KHz (i.e., a time constant ~20 psec). While 
PbS detectors have a time constant ranging from 5 msec (uncooled) to 
0.5 msec (cooled to 77°K). 

In order to meet the Thematic Mapper bandpass requirements an InAs 

detector cooled to 77°K would be required. This would provide a time 

constant Jslpsec and a D* value of around 3 x 10^^. Although this value 

1 

of detectivity is marginal, it could be improved by cold stopping or 
cold filtering the detector. The corresponding noise-equivalent power, 
NEP, IS given by: 

A 

( 6 ) 

using B = 51.2 X 10^ hertz, A^ = 0.0104 cm^, and D* = 3 x 10^^ gives 
NEP = 7.84 X 10“^^ watts. 
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Neglecting the attenuation of the radiation reflected by the 
ground scene as it passes from earth to the sensor and through the 
optxcal system of the sensor, the expected signal in this band from a 
material of 50% spectral reflectance would be 3.72 x lO”^® watts. Hence, 
the corresponding noise equivalent change m reflectance, NEAp, would 
only be w2%. 

9.4 Conclusions 

Evidence for incorporating a band between 1.2 and 1.3 micrometers 
(there are water absorption bands at 1.15 and 1.45 micrometers) is 
inferential and rests upon the idea that having two observations in the 
near infrared would permit a resolution element to be better character- 
ized in terms of a mixture of vegetation and soil. This band was not 
selected because low detector sensitivity combined with the already 
high performance of the TM in agricultural applications made this 
spectral band a 3rd choice. 
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10 

CANDIDATE BAND IN THE 2.2 pm REGION 


10.1 Background 

Of particular interest to geologists is the 2.0-2.6jim region, 
which contains absorption bands caused by electronic processes associ- 
ated with the hydroxyl radical (OH) , sulfate radical (SO^) , carbonate 
radical (CO^) and a few less important absorption bands in rocks and 
minerals (Hunt, G. B., 1977). A cellulose absorption band also occurs in 
this region and is thought to be useful for vegetation mapping among 
diseased or dead vegetation. This explains, why a number of scientists 
from different fields are interested in this wavelength region as a 
candidate for the 7th band. 

For mapping limonitic altered rocks , the current 4-channel Landsat 
MSS sensor is very effective because it responds to the intensity of 
ferric-iron. The current MSS is ineffective for mapping nonlimonitic 
altered rocks. Also, nonhydro thermally altered limonitic rocks (such as 
limonitic tuffs and flows) are commonly indistinguishable from limon- 
itic rocks. By recording multispectral images at longer wavelengths 
than is possible with the Landsat MSS, these limitations might be over- 
come (Abrams, M. J., et al, 1977). 

Of greatest interest is the 2.15-2.25ym region, where nearly all 
altered rock spectra exhibit an intense absorption band (referred to 
here as the 2.2ym dip). Although the physical origins of this dip at 
2.2ym are quite complex (see Appendix A) this spectral feature is largely 
independent of limonlte content. On the other hand this 2.2ym absorp- 
tion dip is absent or weak in the spectra of most unaltered rock. 

Figure 8 shows typical spectral reflectance curves for altered and 
unaltered rocks. It should be noted that the limited solar energy 
available in this region precludes adding a Thematic Mapper band which 
will unequivocally respond to this 2.2ym dip. In addition, care must 
be taken not to pick a spectral band which would allov? the 2.35ym dip 
in carbonates to be confused with the 2.2pm dip of altered rocks. 
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(b) 


FIGURE 8. SPECTRAL REFLECTANCE CURVES OF THREE ALTERED (a) AND 
THREE UNALTERED (b) ROCKS FROM THE GOLDFIELD, NEVADA 
MINING DISTRICT (after L. Rowan et al, 1977) 
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For most altered roclcs the reflectance maximum is at about 1.7]jm; 
the reflectance level at 2.2ym is relatively low compared to Thematic 
Mapper channel 5 tl. 55-1. 75ym) , whereas reflectances at these two wave- 

S 

lengths are similar for most unaltered rocks. In order to determine 
the usefulness of these spectral differences for discriminating non- 
limonitic altered rocks from unaltered rocks not only the absolute 
reflectances but also the ratios of the spectral bands can be analyzed. 

For that reason it was ascertained by R. Vincent (1977), that it 
would be theoretically possible to separate altered from unaltered 
rocks with currently planned Landsat-D bands. He makes this statement 
based on an article published by L. Rowan et al, 1977 where rock and 
mineral samples are discussed that were collected in the Goldfield 
Nevada Mining District. The spectral reflectance curves shown in 
Figure 8 were used by Rowan to illustrate the usefulness of a 2.2ym 
spectral band in discriminating zones of hydro thermally altered rocks 
from unaltered rocks. However, as was subsequently pointed out by 
R. Vincent for these specific spec era it is possible to separate the 
three altered from the three unaltered rocks using ratios of currently 
planned Landsat-D spectral channels. 

In response to this suggestion, recent work at NASA/JPL using a 
large number of representative rock spectra has shown that the current 
Landsat-D Thematic Mapper bands are limited in their ability to dis- 
criminate altered from unaltered rocks (J. M. Soha, private communica- 
tion, September 26, 1977). 

It is well kno\im that natural vegetation can significantly mask 
and alter the spectral response of the ground, and so it is likely, that 
vegetative cover will always cause some problems in this wavelength 
region for geologic remote sensing. However, an additional band in this 
region would be quite useful for geology and would also help improve 
the discrimination among various classes of vegetation with different 
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degrees of vigor. The usefulness of a spectral band in this wavelength 
region might be reduced by several considerations from physical and 
engineering standpoints : 

-2 -1 

1. The solar irradiance at sea level (W m pm ) is diminished 
by a factor of about 25, compared to the maximum in the green 
wavelength region. 

2. In addition, the solar radiation and the emitted radiation 
from the target Itself are approximately equal near the 
"cross-over" in the 2. 0-2. 6pm region. It is expected that 
this may result in a lower signal- to-noise ratio than could 
be achieved in a shorter wavelength region. 

3. Since the addition of a 7^^ channel in this wavelength region 
would require an independent detector dewar package, this 
would place a significantly increased demand on the sensor 
support system. 

10.2 Quantif ication 

10.2.1 Introduction 

In order to provide a quantitative measure of the usefulness of 
a band centered near 2.2pm a linear discriminant analysis was performed 
on a large number of field and laboratory spectra. The majority of 
these spectra were provided to us by NASA/JPL through the good offices 
of our technical monitor. The calibrated field spectra, which were 
collected by NASA/JPL with a Portabde Field Reflectance Spectrometer 
(PFRS) during the years 1974-1977 at a wide variety of sites in the 
Western United States, formed an excellent technical basis for this 
study. These, yet unpublished, spectra were grouped as shown in Fig- 
ure 9 . 
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Altered Rocks Mo . of Samples 

Limonitic Altered - <5% iron 10 

Hematitic altered - >5% iron 87 

Non-Iron Bearing Altered - bleached 13 

Unaltered Rocks 

Rhyodacite 15 

Andesite 29 

Rhyodacite Ash Flow Tuff 10 

Rhyolite Ash Flow Tuff 11 

Anorthosite 20 

Granite 18 

Red Sandstone 8 

Red-Pink Shale-Mudstone 5 

Basalt 14 

Coaldale Red Tuff 20 

Coaldale Light Tuff 62 

Argillite 12 

Carbonates 29 

Total Number of Samples 363 


FIGURE 9. GROUPING OF ROCK SPECTRA FROM NASA/JPL 


To these 363 rock spectra from NASA/JPL we added additional spectra 
from the NASA Earth Resources Spectral Information System (ERSIS) to 
assure that from a general geologic point of view all important class 
types were well represented. Figure 10 indicates the breakdom of the 
78 additional spectra from ERSIS. 
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Clay Soils 2 

Loam Soils 13 

Sandy Soils and Sand 7 

Carbonates (Minerals) 26 

Carbonates (Rocks) 3 

Sulphates 3 

Acidic Silicate Rocks 6 

Intermediate Silicate Rocks 5 

Basic and Ultrabasic Silicate Rocks 6 

Healthy Vegetation (Corn, Sorghum, Grass, 

Clover, Pine) 5 

Dead Vegetation (Corn, Sorghum) 2 

Total Number of Samples 78 


FIGURE 10. GROUPING OF ADDITIONAL SPECTRA FROM THE 
EARTH RESOURCES SPECTRAL INFORMATION SYSTEM 

Performing a linear discriminant analysis on the average inband 
reflectance values for the 363 samples from NASA/JPL we were able to 
discriminate the altered from unaltered rocks with an average accuracy 
of 78% using the current 5 TM-Bands. By including the additional band 
which was recommended by NASA/JPL (2.08-2. 35ym) , we obtained an average 
accuracy of 96% correct classification. 

10.2.2 Candidates 

As mentioned in Section 4, the different opinions of several 
scientists concerning the best spectral bandwidth, required more detailed 
examination. Different considerations had to be taken into account: 
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- Well positioned with regard to the different absorption 
bands. Candidate bands were only selected from the region 
between 2.0 and 2.5ym because of the atmospheric water 
absorption band at 1.9pm and the combined water and carbon 
dioxide absorption bands which make the atmosphere opaque 
in the region just above 2.5ym (see Figure 3). 

- Bandwidth as narrow as possible relative to sensor noise. 

The different candidate spectral bands (referred to here as Channels 6-13) 
which were investigated are shown in Figure 11. 


Channel Number Spectral Interval 


Channel 

6 

2.08 - 

2.35 

ym 

Channel 

7 

2.15 - 

2.35 

ym 

Channel 

8 

2.10 - 

2.35 

ym 

Channel 

9 

2.05 - 

2.35 

ym 

Channel 

10 

2.10 - 

2.40 

ym 

Channel 

11 

2.05 - 

2.45 

ym 

Channel 

12 

2.00 - 

2.40 

ym 

Channel 

13 

2.00 - 

2.50 

ym 

FIGURE 11. 
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10.2.3 Transfer to a Spaceborne Sensor 

To determine which of these candxdate spectral bands would provide 
the best performance over the wide variety of conditions of observation 
which will be encountered by the Thematic Mapper, as it views the Earth 
from arctic to the tropics, computer simulations were required. First 
the inband radiances at the top of the atmosphere were computed for all 
441 samples under a given condition of observation. Then these inband 
radiances were quantized in a manner which simulated the presence of sen- 
sor noise, before a linear discriminant analysis was carried out to pre- 
dict the performance of a particular spectral configuration under this 
condition of observation. This entire procedure was repeated for several 
different atmospheric conditions . 

10.2.3.1 Atmospheres/Conditions of Observation 

For our computations we chose to use a set of seven different 
atmospheres /conditions of observation which are representative of the 
range of situations the Thematic, Mapper will encounter in supporting a 
broad community of users. These conditions are shown in Figure 12. 

Cases A-F range from the excellent viewing conditions obtained when 

_ , r 

viewing the nadir at high noon through a tropical mid-day atmosphere with 
a visual range of 23 kilometers ,(case E) to the limited viewing condi- 
tions obtained when viewing the Earth at a high latitude during the win- 
ter time through a temperate mid-day atmosphere with a visual range of 
only 13 kilometers X'jhile the sun is 67.5° from zenith (case C) . Case G 
was carried along only as an extreme. It corresponds to case C with the 
exception that visual range has been reduced from moderate haze (13 kilo- 
meters) to severe haze (5 kilometers). 
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Aerosol 

Distribution 

Visual 

Range 

Sensor Viewing 
Angle 

Solar Zenith 
Angle 

Case A 

Tropical 

13 

km 

0° 

45° 

Case B 

Tropical 

23 

km 

0° 

45° 

Case C 

Temperate (winter) 

13 

km 

0° 

67.5° 

Case D 

Temperate (winter) 

23 

km 

0° 

67.5° 

Case E 

Tropical 

23 

km 

0° 

0° 

Case F 

Arctic (summer) 

23 

kra 

0° 

45° 

Case G 

Temperate (winter) 

5 

km 

0° 

67.5° 


FIGURE 12. CONDITIONS OF OBSERVATION USED IN PREDICTING 
THEMATIC MAPPER PERFORMANCE 

For the reflective region covered by TM bands 1-5 the radiative 
transfer- model can be exercised once for each condition of observation 
to produce two functions, and t^(X) which are then used to 

transform ground based values of spectral reflectance, p(A), for a 
sample into the corresponding spectral radiance, L(A), at the top of 
the atmosphere. 

After we noted from our computations that the radiant energy emitted 
by a ground material at 300 °K contributed less than 1% of the total radi- 
ance at the top of the atmosphere (consider Figure 3), we employed this 
same technique in the 2. 0-2. 5pm region. 

For a given atmosphere/condition of observation, integration of L(A) 
for each sample across the spectral intervals corresponding to the current 
TM bands and the candidate bands shotm. in Figure 11 produced a thirteen- 
dimensional observation vector. 
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10.2.3.2 Predicted Performance Without Sensor Noise 

Figure 13 and Figure 14 show the predicted performance of the 
Thematic Mapper in discriminating hydro thermally altered rocks from 
unaltered rocks for each of the candidate bands investigated under each 
atmosphere/condition of observation with no simulation of noise or quan- 
tization effects. These results were obtained by applying a linear 
discriminant analysis directly on the thirteen-dimensional inband ra- 
diance vectors of Section 2.3.1. Column 1 gives the baseline perform- 
ance obtained using the five reflective bands of the current Landsat-D 
TM. The second through ninth columns give the percent correct classi- 
fication which can be expected under this same range of observing condi- 
tions by inclusion of candidate spectral band 6 through candidate spec- 
tral band 13, respectively. The most significant item in Figure 13 
(363 samples from NASA/JPL) is that the predicted performance gain is 
essentially the same for all candidate bands and does not even vary 
with the atmospheric conditions. The reason for this result is that, 
in the absence of sensor noise or quantization, the effect of the ra- 
diative transfer model for the reflective region is to transform the 
inband reflectance of each sample into a corresponding radiance value 
through essentially an affine transformation whose coefficients are deter- 
mined by the specific condition of observation (i.e.. Case A, Case B, 
etc.). Hence, the double precision linear discriminant analysis which 
was used rescaled the variables involved and achieved the same degree 
of separability in all cases. When the 78 additional spectra from 
ERSIS are included (Figure 14) there was an overall decrease in 
correct classification with small differences from atmosphere to at- 
mosphere but a similar result was obtained. 

To overcome this limitation the effect of sensor noise had to be 
taken into account . 
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1-5 

+6 

+7 

+8 

+9 

+10 

+11 
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A 
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96.2 
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94.7 
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B 

79.4 

95.1 

96.2 

95.1 

95.1 

95.1 

95.1 

94.7 

94.7 

C 

79.4 

95.1 

96.2 

95.1 

95.1 

95.1 

95.1 

94.3 

94.1 

D 

79.4 

95.1 

96.2 

95.1 

95.1 

95.1 

95.1 

94.3 

94.3 

E 

79.4 

95.1 

96.2 

95.1 

95.1 

95.1 

95.1 

94.5 

94.3 

F 

1 

79.4 

95.1 

96.2 

95.1 

95.1 

95.1 

95.1 

94.3 

94.3 

Mean 

79.4% 

95.1 

96.2 

95.1 

95.1 

95.1 

95.1 

94.5 

94.4% 

Std. Dev. 

o 

o 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.18 

0.24% 


FIGURE 13. PREDICTED PERFORtlANCE OF THE THEMTIC MAPPER IN DISCRIMINATING 
HYDROTHEEMALLY ALTERED ROCKS FROM UNALTERED ROCKS 
(Percent Correct Classification) 

(No Sensor Noise — 363 Samples) 
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Illumination 

Spectral 

Channels 

1-5 
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+8 

+9 

+10 

+11 

+12 

+13 

A 

76. IX 

93.1 

93.4 

93.2 
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92.6 

92.9 
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B 

76.4 

93.1 

93.4 

93.2 

92.8 

92.4 

92.6 
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C 
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92.8 
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E 
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94.0 
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93.1 

93.1 

F 
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92.9 

93.4 

93.0 

92,5 
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92.3 
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Mean 

77.40% 

93.17 

1 

93.73 

93.43 

93.07 

92.95 

92.72 

92.87 

92.88% 

Std, Dev. 

1.31% 

0.20 

0.64 

0.37 

0.36 

0.51 

0.39 

0.12 

0.12% 


FIGURE 14. PREDICTED PERFORMANCE OF THE THEMATIC MAPPER IN DISCRIMINATING 
HYDROTHERIIALLY ALTERED ROCKS FROM UNALTERED ROCKS 
(Percent Correct Classification) 

(No Sensor Noxse — 441 Samples) 
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10.2.3.3 Sensor Noise 

As is well known, the performance of a sensor can be defined rel- 
ative to parameters such as IFOV, spectral resolution and signal- to- 
noise (S/N) ratio. The noise in a sensor may arise from the radiation, 
the current induced by the incoming radiation, the detector itself or 
the detector electronics but the predominant source is dependent upon 
the detector and its operating conditions . 

For the reflective region the noise equivalent change in radiance, 
NEAL, depends upon the signal level. This dependence is well represented 
by a function of the form: 

NEAL = /a + gS (7) 

where S is the signal level while a and 6 are constants which characterize 
respectively the dark current of the detector and the shot noise of the 
detector due to both the statistical nature of the local photons and the 
response of the detector to the incoming electromagnetic radiation. 

Using the sensor performance figures (S/N ratios at the maximum 
and minimum radiance levels) given in the Landsat-D Thematic Mapper 
sensor specifications (GSFC, 1977) we calculated values for a and 6 for 
each of the existing spectral bands. Since an InSb detector was the 
best choice for a 2.2)Jm centered spectral band we used the values of a 
and 3 from the existing 1.55-1.75ym band which uses this same detector 
type. Figure 15 shows the values of a and g we obtained and the 
corresponding values of NEAL at the expected minimum and maximum radiance 
levels. As can be seen, for all of these bands we are operating near 
the knee of the curve between a dark current limited sensor and a shot 
noise dominated sensor. For our computer simulation of sensor noise 

we used quantization levels which corresponded* to a one sigma NEAL 

-2 -1 -2 -1 
value of 10 ywatts cm sr for channels 1-4 and 7 ywatts cm sr 

NEAL for channels 5-13. 


*See Appendix C 
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FIGURE 15. PARAMETERS DESCRIBING THE NOISE EQUIVALENT RADIANCE LEVELS 
OF THE LANDSAT-D THEMATIC MAPPER CHANNELS 
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10.2.3.4 Predicted Performance with Sensor Noise Included 

Figure, 16 shows the predicted performance of the Thematxc Mapper 
in discriminating hydrothermally altered rocks from unaltered rocks 
under a variety of conditions of observation with the inclusion of la 
sensor noise. These results are based upon a linear discriminant anal- 
ysis of the 441 samples obtained by simulating the effect of sensor 
noise after applying the radiative transfer model to the combined spectra 
from both NASA/JPL and the ERSIS. 

For a given atmosphere/condition of observation the addition of 
a sixth reflective band gives an increase in the percent correct classi- 
fication, with candidate spectral bands 8, 9, 10, and 11 giving geneially 
lower performance than spectral bands 6 (2. 08-2. 35ym) , candidate 7 
(2.15-2.35 pm)j or candidate 13 (2.00-2.50 pm). 

To understand' better why four of the seven candidate bands per- 
formed significantly poorer than the remaining three candidate bands 
we determined which spectra were being misclassif ied . Inspection re- 
vealed that the majority of this difference in performance was due to 
a large proportion of the Non- Iron bearing altered rocks, the Coaldale 
Red Tuff, and the Coaldale Light Tuff being consistently misclassif ied 
when one of _these four spectral bands were used. 

Because the goal was to select a spectral band which gave the best 
possible performance over a wide range of conditions of observation we 
used the mean classification accuracy over all six atmosphere/ illumination 
conditions as our figure of merit. In order to ascertain the dependence 
of these results on the noise level which was modeled, we repeated all 
the computations twice more: 

- Noise Level - a/ 3 

- Noise Level = 3a 

The results of these simulations are given in Figure 17 . 
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FIGURE 16. PREDICTED PERFOPvMANCE OF THE THEilATIC llAPPER IN DISCRIMINATING 
HYDROTHERMALLY ALTERED ROCKS FROM UNALTERED ROCKS 
(Percent Correct .Classification) 

(lo Noise Level — 441 Samples) 
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FIGURE 17. SENSITIVITY OF PREDICTED THEMATIC MAPPER PERFORiMANCE TO SYSTEM NOISE; 
HYDROTHERMALLY ALTERED ROCKS FROM UNALTERED ROCKS 
(Percent Correct Classxf ication) 

(441 Samples) 
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Our xnterpretation of these results is that in the presence of 
small noise levels (jsa/3) the best performance is obtained from the 
narrowest spectral band (Candidate #7, 2.15-2.35 ]im) . In the presence 
of high noise levels ( w3(J) the best performance is achieved by using 
the broadest spectral band (Candidate #13, 2.00-2.50 ym) . '.For a system 
whose noise levels conform to the Landsat-D sensor specification (ssla) 
the 2.08-2.35 ym band is optimum followed closely by the 2.00-2.50 ym 
band. Since the primary physical phenomena being exploited is the 
2.2 ym dip in the reflection spectra of hydro thermally altered rocks, 
we opted for the first of these two candidates on the grounds that all 
other things being equal, a narrow spectral band will provide results 
which are more unequivocal. 

10.2.3.5 Additional Investigations 

One of the technical concerns about a 2.2 ym centered band was the 
influence of the carbonates, i.e., would carbonates become confused 
with hydro thermally altered rocks. To answer this question we examined 
the performance on a special data set of 168 samples (110 hydro thermally 
altered samples, 29 carbonates samples from NASA/JPL and 29 carbonate 
rocks and minerals from ERSIS). Figure 18 shows the predicted per- 
formance of the TM in discriminating between altered rocks and carbonates 
in the presence of Icf noise. 

As can be seen from this table both the 2.08-2.35 ym band and the 
2.00-2.50 ym band perform quite well. Considering that the carbonate 

absorption dip around 2.35 ym is somewhat masked by the rapid decrease 
in solar irradiance with increasing wavelength, it is not surprising 
that the broader band is slightly better and more consistent. 
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FIGURE 18. PREDICTED PERFOEMANCE OF THE THEMTIC MAPPER IN DISCRIMINATING 
HYDROTHERMALLY ALTERED ROCKS FROM CARBONATES 
(la Noise Level — 168 Samples) 
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We also checked the performance of a 2.08-2.35 pm band for general 
geologic purposes where not only altered rocks should be discriminated 
from unaltered rocks, but all rock types among themselves. This was done 
by applying a linear discriminant analyses to tiie 15 band ratios obtained 
for 390 samples in the presence of la noise for both the best and worst 
condition of observation. The results are shown in Figure 19. As 
expected, there is a considerable improvement in performance between 
atmosphere C and atmosphere E because ratio processing is most appro- 
priate when atmospheric path radiance is small. It should be noted 
that for general mapping the average class performance we obtained 
ranged between only 54% correct classification and 62% correct classi- 
fication. 

For completeness, we include a correlation matrix for each of the 
original TM bands and the candidate spectral bands considered in this 
section. This matrix (shown in Figure 20) was computed for atmosphere E 
and a la noise level over all 441 sample spectra. 

10 . 3 Conclusions 

Based upon the substantial improvement in mapping zones of hydro- 
thermally altered rocks from unaltered zones over a broad range of 
illumination and atmspheric conditions, a 2.08-2.35 pm spectral band 
having a ground resolution of 30 meters is recommended for inclusion on 
the Landsat— D Thematic Mapper. Imagery from this spectral band will 
also be helpful for general geologic mapping; as well as for discrim- 
ination among various classes of vegetation with different degrees of 
vigor. 
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Percent Correct 

Classification 

Sample 

Number 

Worst Case 
(Atmosphere C) 

Best Case 
(Atmosphere E) 

Liraonitic Altered 

10 

20 ,0 % 

30.0 % 

Hema title Altered 

87 

70.1 

79.3 

Non- Iron Bearing 

13 

84.6 

84.6 

Rhyodacite 

15 

40-0 

40.0 

Andesite 

29 

44.8 

58.6 

Rhyodacite Ash Flow Tuff 

10 

60.0 

100.0 

Rhyolite Ash Flow Tuff 

11 

54.5 

72.7 

Anorthosite 

20 

35.0 

20.0 

Gram te 

18 

33.3 

50.0 

Red Sandstone 

8 

87.5 

62.5 

Red-Pink Shale-Mudstone 

5 

60.0 

80.0 

Basalt 

14 

64.3 

85.7 

Coaldale Red Tuff 

20 

90.0 

85.0 

Coaldale Light Tuff 
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FIGURE 19. PREDICTED PERFORMANCE OF THE THEMATIC MAPPER 
FOR GENERAL GEOLOGIC MAPPING 
(Ratio Processing, 390 samples. 

Channels 1-5+6, Icr noise) 
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FIGURE 20. CORRELATION JIATRIX FOR THE THEMATIC MAPPER CHANNELS USED IN THIS SECTION 

(Atmosphere E, Iff Noise Level — All 441 Samples) 



O O 

wp 

<o 

c; 

^ s? 



FORMERLY WILLOW RUN LABORATORIES THE UNIVERSITY OF MICHIGAN 



11 

CANDIDATE BAND IN THE THERMAL IR REGION 
11.1 Introduction, Background 

It IS assumed that multispectral thermal infrared scanners with 
more than one channel in the 8-14 pm wavelength region are xmportant 
for geological remote sensing because the stretching and rotational 
processes of SiO^ tetrahedra in silicate rocks and minerals produce 
emittance minima in this wavelength region which occur at different 
wavelengths, depending on the mineral composition (R. Lyon 1964). The 
reasons why emittance minima of the silicate minerals occur at different 
wavelengths are very complex (see R. Vincent 1975a) and cannot be dis- 
cussed here. By contrast, vegetation shows nearly a constant emittance 
in that spectral region and most nonsilicate rocks and minerals exhibit 
a quite different emittance minima from silicates. For instance, 
carbonate rocks (limestone, dolomite) have only one rather sharp band 
near 11.3 pm. 

R. Vincent (1973b) has theoretically examined this systematic shift 
of emittance minima to longer wavelength with decreased quartz content 
and increased iron and magnesium content. Working with low-altitude air- 
craft multispectral scanner data he also demonstrated (R. Vincent, 1974 
that the optimum filter combination for the purpose of discriminating 
among silicate rock types on the basis of % 8 IO 2 and % (S 1 O 2 - Al^O^) was 
8.1-10.1 pm for the shorter and 9.2-11.2 pm for the longer wavelength 
filter. This division of the 8-12 pm region into two thermal channels 
is supposed to enhance emittance variations in exposed silicate rock 
outcrops, while suppressing temperature variations across the scene. 

From the many results achieved by researchers seeking to exploit 
the supposed two -channel- thermal IR capabilities for geologic applica- 
tions only a few will be enumerated: 

- using a two-channel thermal IR ratio a rhyolitic tuff was 

discriminated from alluvium (Vincent 1975a). (This raises the 
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question whether the unaltered pink rhyolitic tuff mentioned 
by Rowan et al 1977 could have been separated from altered rocks 
in the Goldfield area using a thermal IR ratio processing. 

Because of the great difference in quartz content between 
these rock types, it is likely, that the tuff and the montmoril- 
lonitic-argillic rocks could have been discriminated with a two- 
channel IR ratio image) , 

- in general a two-channel IR capability would assist in dis- 
criminating alteration zones that are typified by silicification 
in background rocks of low quartz content and kaolinization from 
background rocks of high quartz content . The ability to dis- 
criminate .silicates from nonsilicates would be useful in tungsten 
exploration. 

- two-channel IR ratio processing would permit discrimination 
between sandstones and limestones or other carbonate rocks as 
well as quartz sand from calcite sand in beach reconnaissance. 

- the mapping of ultrabasic silicates related to diamond-producing 
kimberlites such as ophiolites should be readily performed with 
a two-channel IR scanner. 

- even sulphates and phosphates, which also have an emittance mini- 
mum in this wavelength region and .compete with silicates, could 
be discriminated by ratio processing. 

In addition to these geological applications the following 
reasons for supporting a second thermal IR channel have also been 
suggested: 

- two thermal bands would provide a backup band .for thermal 
information so that if one channel malfunctions, the other 
could be used for such measurements,, even though the two-channel 
ratio method could no longer be performed ._ 

“ by exploiting the differences in emission it should be possible 
to discriminate dead and dry vegetation. 
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- thermal band cooling requirements would probably represent less 
of an increase in the demand upon the sensor support system 
than that required for a 2 . 2iim centered band provided the 2nd 
themal sensor was placed in the same dewar. 

11.2 Quantization 

11.2.1 Introduction 

In order to provide a quantitative measure of the usefulness of a 
second thermal IR channel for geological application, we began by making 
use of the ERSIS data. It was not as easy as one might assume to find 
measured rock spectra in the 8-14iJm wavelength region which are mineral- 
ogical more or less unequivocal. Unfortunately, we were not able to 
find as many samples of rock spectra as we would have preferred to 
support a quantitative study such as this. Also, no comparable study 
of discrimination between altered and unaltered rocks could be conducted 
using the thermal region due to lack of appropriate emission spectra. 

It would have been best to use the same rock samples or at least the 
same kind of rocks for both the reflective and emissive portion of this 
study, but this was not possible. Figure 21 shows the different classes 
which could be established, and the number of individual spectra which 
each class contained : 

No. of 


Class Spectra 

Acidic Silicate Rocks 47 

Intermediate Silicate Rocks 16 

Basic and Ultrabasic Silicate Rocks 20 
Carbonate Rocks 13 

Loams 8 

Silicate Mineral & Silicate Sand 4 

Hydro thermally Altered Rocks 3 

Total Number of Samples 111 


FIGURE 21. GROUPING OP ROCK SPECTRA IN THE THERMAL REGION 
FROM THE EARTH RESOURCES SPECTRAL INFORtl&TION SYSTEM 
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ri.2.2' Transfer to Space 

As explained earlier the radiative transfer model also took into 
account solar illumination, atmospheric absorption and attenuation in 
the thermal IR region. For economy we limited ourselves to only one 
type of atmosphere, namely: 

Case B: Tropical atmosphere, visual range = 23 km, 

sensor view angle = 0°, solar zenith angles = 45°. 

Because of the ozone absorption band in the 9.4-9. 6 ym region, it 
was suggested (Vincent, 1977) that the spectral bands be located at 
8.2-9. 2 ym and 9.8-11'. 8 ym instead of the overlapping bands from 
8.1-10.1 ym and 9.2-11.2 ym which were used in his 1974 study. 

11.2.3 Signal/Noise Considerations 

This investigation also required that attention be given to the 
signal/noise characteristics of the Landsat-D Thematic Mapper. We used 
the following value: 

NEAL = 6.40 yW sr”^cm"^ which corresponds to NEAT of 0.5°K at T=300°K 
and simulated sensor noise by the procedure presented in Appendix C. 

11.2.4 Classification Results 

We defined the following three spectral bands for use in the ther- 
mal IR study with 111 samples under atmosphere/ condition of observation B; 


Channel X = 10.4-12.4 ym 
Channel Y = 8 . 2-9 . 2 ym 
Channel Z = 9.8-11.8 ym 

The performance of the Ttf for general geologic mapping using these thermal 
channels alone is shown in Figure 22 for the case of no sensor noise, and 
in Figure 23 for the case of lo sensor noise. In the configuration where 
tx^o thermal channels are available, we also show results using ratio 
processing. 
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No Sensor Noise — T 

=300 '■K 


X 

Y & Z 

Y/Z 

Acidic Silicate Rocks 



48.9% 

51. 

.1% 

25.5% 

Intermediate Silicate 

Rocks 


12.5 

37. 

,5 

31.3 

Basic and Uitrabasic 

Silicate 

Rocks 

10.0 

15. 

,0 

10.0 

Carbonate Rocks 



15.4 

46. 

.2 

61.5 

Loams 



100.0 

100. 

,0 

37.5 

Silicate Mineral and 

Silicate 

Sand 

50.0 

75. 

.0 

75.0 

Hydrothermally Altered Rocks 


33.3% 

66. 

.7% 

66.7% 


FIGURE 22. PERFORMANCE OF THEMATIC MAPPER THERMAL CHANNELS 
FOR GENERAL GEOLOGIC MAPPING (NO NOISE) 


la .Sensor Noise — T=300°K 


X 

Y & Z 

Y/Z 

Acidic Silicate Rocks 


57.4% 

53.3% 

21.3% 

Intermediate Silicate Rocks 


6.3 

43.8 

12.5 

Basic and Uitrabasic Silicate 

Rocks 

— 

15.0 

55.0 

Carbonate Rocks 


46.2 

46.2 

30.8 

Loams 


100.0 

100.0 

37.5 

Silicate Mineral and Silicate 

Sand 

50.0 

75.0 

75.0 

Hydrothermally Altered Rocks 


33.3% 

66.7% 

66.7% 


Mean 

42% 

57% 

43% 


FIGURE 23. PERFORMANCE OF THEMATIC MAPPER THERMAL CHANNELS 
FOR GENERAL GEOLOGIC liAPPING (la SENSOR NOISE LEVEL) 
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Both' of” these simulations were conducted using the same temperature 
(300 °K) for 'all the samples. We note that: 

(a) the classification accuracy is far below the authors’ 
expectation 

(b) the inclusion of a second thermal IR channel did not 
provide a great increase in classification accuracies. 

(c) performance of the ratio processing is disappointingly 
low in the case of la sensor noise. 

(d) the performance which can be expected from a single 
thermal channel is very limited. 

Because the use of a single temperature for all the samples is unre- 
alistic even for observations at 9:30 AM local time, we repeated the 
simulation with the temperature of each sample drawn from a gaussian 
distribution. The idea behind this was that for a given scene, there 
would be no dramatic differences in temperatures from class to- class but 
the samples would range in temperature according to different illumina- 
tion conditions . 

In accordance with this reasoning we drew 111 temperatures from a 
normal distribution which had a mean value of 300°K and a standard devi- 
ation of 5“K. The radiative transfer model was then exercised using one 
of these temperatures for each sample. The predicted performance of the 
TM thermal channels under these conditions is shown in Figure 24. 
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la Sensor Noise 


X 

Y & Z 

Y/Z 

Acidic Silicate Rocks 


21.3% 

40.8% 

44.7% 

Intermediate Silicate Rocks 


25.0 

'56.3 

37.5 

Basic and Ultrabasic Silicate 

Rocks 

20.0 

35.0 

40.0 

Carbonate Rocks 


— 

30.8 

7.7 

Loams 


87.5 

75.0 

12.5 

Silicate Mineral and Silicate 

Sand 

0 

75.0 

75.0 

Hydrothermally Altered Rocks 


33.3% 

33.3% 

66.7% 

Mean 


50% 

41% 


*Mean sample temperature = 300°K with a standard deviation of 

FIGURE 24. PREDICTED PERF0RI4MCE OF THEtlATIC MAPPER 
THERMAL CHANNELS FOR GENERAL GEOLOGIC MAPPING 
(Random* In-scene Temperatures) 

Our interpretation of these results leads to the following state- 
ments : 

(a) Compared to the performance with two thermal channels, a 
single thermal channel provides poor results, 

(b) For general geologic mapping into these seven classes, two 
thermal channels gave classification accuracies close to the 
accuracies achieved by the TM with six reflective channels on 
the eighteen classes used in Section 10. 

(c) The somewhat lower classification accuracies achieved using 
ratio processing is indicative of the decrease in classification 
accuracy which should be expected for any preprocessing proce- 
dure which accounts for the scene dependent external effects 

in thermal data. 

(d) Further study needs to be made to identify the most effective 
channels in thermal region for general geologic mapping. 
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11.3 Conclusions 

We emphasize that this thermal study represents no comparable study 
to the one performed in the 2.2 pm regxon due to the limited amount of 
spectra which are available. Incorporation of a 2nd thermal band 
(8. 2-9. 2 micrometers) would improve classification accuracies for some 
general rock types over that obtainable with a single thermal band. 
Although the present work does not support the contention that accurate 
identification of rocks according to their S 1 O 2 content would be possible 
with two thermal bands, a single thermal band is inadequate to separate 
variations due to conditions of observation from changes in ground tem- 
peratures and emissivities . The authors do not agree with the remark 
mentioned earlier, by a remote sensing scientist that the implementation 
of a 2nd thermal IR band would be premature and should wait until expe- 
rience was gained with Landsat-C data. We fear that although the 
thermal region contains substantial new information, the difficulties 
encountered in trying to use a single channel to account for the three 
competing physical phenomena will Introduce more .problems than it solves. 
In addition to these difficulties the current 120 meter ground resolution 
of the TM thermal band will require great care in the way these thermal 
data are used in conjunction with the reflective bands during machine 
processing. Also further studies are recommended to penetrate the 
relationship between the thermal and reflective region as well as the 
potential of the thermal region itself for general geologic mapping. 

Because the performance of an additional thermal band for general 
geologic mapping is on the same order as that obtainable with a 2.2 pm 
band, this candidate spectral band is considered to be a 2nd choice. 
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12 

COMPARISON WITH CURRENT LANDSAT MSS 

In order to determine the usefulness of adding different spectral 
bands to improve the capability given by the present 4-channel Landsat 
MSS, we ran a study using the same rock spectra collected by NASA/JPL 
as well as those from the ERSIS. We took into consideration the sensor 
specification for the current Landsat-II channels and used a linear 
discriminant analysis mth the best case of atmospheric condition 
(atmosphere E) and a lo signal/noise quantization to distinguish altered 
from unaltered rocks. The following accuracies illustrate the perforra- 
ance which can be expected based upon the classification for 440 of our 
samples: 


Landsat-II MSS 
4 Channel 


Current TM 
Channel 1-5 


Recommended TM 
Channel 1-5 (+6 ) 


72.5% 


79.5% 


93.8% 


Statements made to us by the GEOSAT committee, confirm that the 
ability to map zones of altered rocks is very important economically. 
Hence, the dramatic improvement in classification accuracy for zones of 
hydro thermally altered rocks gained by incorporating a 2.08-2.35 ym 
spectral band on the Thematic Mapper justifies its inclusion. 
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APPENDIX A 

SPECTRAL FEATURES OF SILICATE AND CARBONATE ROCKS 
IN THE 2 TO 3 ym REGION 

by 

George Lindquist 


We have done some investigating to try to extract the phenomenology 
of the spectral features seen around 2.3 ym. I believe the explanation 
given in three papers by Hunt and Salisbury (Refs. 1, 2 , 3) is basically 
correct, but some of the identifications given there are somewhat mislead- 
ing. Perhaps a restatement of what is meant is in order, which may serve 
to put the overall situation in a more physically complete perspective. 

I. Silicates 


Consider first the water molecule and OH hyrdoxyl complex. 
The spectrum of the complexes are well known. The water 
molecule looks like u 


0 

and has 3 modes of vibration: a syranietric_ stretch at 3654 cm ^ 
(2.74 pm), a bending vibration at 1595 cm~^ (6.27 pm), and an 
assyrometric stretch at 3755 cm“^ C2.66 pm). For this bent mole- 
cule, all these vibrational modes lead to infrared bands. The 
OH complex is a simple diatomic and can hence have only one 
vibrational stretch mode; it is at 3652 cm“^ (2.74 pm), right on 
top of the symmetric stretch vibration of H 2 O. Because of this 
coincidence, it is very difficult to ever see the vibrational 
fundamental band of OH. 



I>Jhen these molecules are bonded to other atoms or complexes , 
as in a crystal, their vibrational motion is hindered somewhat 
(but usually not severely) . This leads to a decrease in their 
vibrational freq^uencies by a few percent which correspondingly 
moves their fundamental bands to longer wavelengths by the same 
percentage. Thus, in liquids and. solids we expect to see the same 
bands from H 2 O and OH as we see in the gaseous state, except that, 
where the orientation of the molecules is fixed, like in an or- 
dered crystal, we expect the bands to be very much narrower than 
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in gases because no rotational structure will be present. Note 
that this is not the situation in liquids however, where the 
molecules are desordered and can interact strongly with each other 
in a variety of orientations. This leads to much broader bands 
in liquids than in either gases or crystals. 

The bands observed in Ref. 1 at 1.4 and 1.9 ym are due to 
the H 2 O molecules and OH complexes directly. The 1.4 pm band is 
the first overtone of OH and/or a combination band of H 2 O 
(corresponding to a simultaneous fransition in both the symmetric 
stretch mode and 'the asymmetric stretch mode). The existence 
of these combination bands, like the existence of the overtone 
bands comes about because of non-linear restoring faces in the 
vibrations and non-linear interactions between the vibrations -- 
they cannot be explained by the simple springs-and-weights analogy 
These combination and overtone bands occur in all molecules . The 
1.9 pm band is a combination band involving a simultaneous tran- 
sition in the bending mode and the asymmetric stretch of H 2 .O. 

These bands (and a much weaker one at 1.1 pm), as well as the 
fundamental obscured by atmospheric water vapor at 2.7 pm, are 
the only bands due solely to water and OH in. the spectral range 
from 1 to 3 pm. ""The presence of a 1.9 "pm band is sure indication 
of the presence of the H 2 O molecule as opposed to the 01! complex 

Whenever H 2 O or OH are attached to other structures, e.g. 
when an OH replaces an 0 in a silicate complex, or when H 2 O 
becomes loosely bonded to a molecule in a crystal; other xnfrared 
bands become possible. These bands would be combination bands 
with one or more of the vibrations associated with the structure . 
to which the OH' or II 2 O is, bonded. These other vibrations could 
be possibly the fundamental of' the silicate ion, ox possibly 
vibrations of the lattice structure in which the molecules are 
arranged or vibrations associated X7ith other degrees of freedom 
in the crystal. The frequency j and hence wavelength', at x^hich 
these combination bands appear must correspond to the sum of the 
frequencies of the modes involved in the transition'^’''. 

The persistent hands seen in powdered silicate in Reference 
1 at 2.1 to 2.3 pm appear only when OH complexes replace one or 
more oxygens, at least according to the chemical descriptions 
included with Ref. 1. This indicates that the bands must be 


'This statement assumes that the energy levels, of the 
vibrational modes are equally spaced which is approximately 
true for the low levels involved at room temperature. 
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combination bands of OH and some other vibration in the molecules. 
From the spectral position of the bands and the fundamental 
frequency of OH, the other moles involved must be between 9 and 
14 ym, viz 


2.1 pm 

OH Fundamental 


4762 cm2 
3652 cm 


Other mode 


1110 cm ^ = 9 pm 


or 


2.3 pm 

OH Fundamental 
Other mode 


4348 cm” ^ 
3652 cm“^ 


695 cm ^ =14.3 pm 


These bands must be combinations between the OH fundamental and 
some vibrational mode whose fundamental is between 9 and 14 pm. 
This frequency is too high for lattice vibrations but lies 
right where one would expect bands of the silicate complex and 
Al-O bands. The strong quartz (Si02) band is right at 9 pm 
(see Figure 1). The Si-0 band, considered as a free complex 
is at 1247 cm” ^ (8,05 pm) . A fully-bonded Si 04 silicate ion has 
a single Si-O stretching fundamental at 555 cm” ^ CH&f* 5 and 6) 

(18 pm). Mineral silicates have been shown to have a variety 
of crystal structures (Ref. 7) , in which oxygen atoms are paired 
with from one to four silicon atoms. Despite the variety of 
crystal structures we might assume that Si-0 stretching funda- 
mentals will all fall between 8.05 pm and 18 pm. On this basis 
the possible combination bands with OH must all fall between 2.04 pm 
and 2.37 pm. This is a relatively narrow spectral range and one 
could postulate that any silicate in which OH replaces one or more 
of the oxygen atoms will have a band in this range. The replace- 
ment of silicon with aluminum and magnesium atoms will modify 
this situation somewhat, but not greatly; both Al and Mg have 
masses similar to Si. 

The exact spectral position and the fine structure of these 
OH- silicate combination bands will be indicative of the parti- 
cular structure of the silicate crystal, including the possible 
substitution of Mg and Al.at the Si sites. 

Hunt (Ref. 3) gives a classification of the bands ob- 
served in several minerals based on these combination bands and 
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and combination bands of, altaminum and,magnesiiam bearing complexes, 
I see no reason to doubt h.is classification at p^resent (see 
Ref, 3,. Figure 8 and pp 5Q8-510X. ' - ^ 

One can therefore make the following statements about 
the signature of silicates; A band between 2 and 2.4 pm in a 
silicate is positive indication of the, presence of OH in the 
silicate. The precise location and structure of the band , 
is indicative of the particular crystal structure. In 'addi- 
tion, in the presence of OH, a b^.nd-will ,also be present at 1.4 
ym. When water molecules are present, bands will be evident at 

I, 4 and 1.9 pm , Very wide bands, such as are present in Figure 

2 , are indicative of liquid water rather than bonded molecules 
These liquid water bands could easily obscure at least the 1.4 pm 
OH band. . . 

II. Carbonates • : 


The carbonates demonstrate features in the 1.5 pm to 2 .'5 pm- 
region which' are characteristic of the carbonate complex (Ref. 

2). The CO^-'-complex has. a symmetric planar ' structure . There 
are six Vibrational modes but- two are doubly- degenerate, meaning - 
that each of the degenerate frequencies arises from two indepen- - 
dent types of motion. . The frequencies and motions of the Car-'- 
bonate complex are shown in Figure 4. Hunt' has described the” 
possible ^combination bands (see Table 1. taken from Ref. 2)'. 

In addition he mentions splitting of these bands due to 
"lifting 'of the' degeneracies"' and/or due to long wavelength - 
vibrational frequencies.. These two, explanations can, be taken, 
as identical; the only.wayc the" degeneracy could be lost would be 
by the presence of asymmetric bonds to the .C 03 --complex which 
would indeed manifest themselves in. the' spectra by combination 
bands with the' vibrations of the CO 3 itself. Thus Hunt's 
explanations of the carbonate ' spectra seem quite reasonable. 

As in silicates details of the carbonate ' spectra will be indica- 
tive of the particular crystals in which the carbonate complexes " 
are involved.' 

III. Signature Implications 

dt ' 

Thus it does appear that both OH-bearing silicates "and 
carbonates will have persistent bands in the 2-2.5 pm region. 
Furthermore the details of these bands will be dependent' on 
the particular crystal being observed. Unless it is possible 
to do a very careful analysis of the spectral details I would 
expect that it would be difficult to reliably, distinguish one 
from the other in remotely gathered images. 
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KA0LIMITE9 MESA ALTA. N M . 



IN MICRONS 


KAOLINITE 2I6B MACOf^, GA. 



FIGURE 1. SPECTRA OF KAOLINITE POWDERS (TTTO SA24PLES) NOTE 
SIMILARITY OF SPECTRAL DETAILS IN 1.4 and 2 , 1 ym 
FEATURES. These details are the manifestation of 
the crystal structure in the OH and OH-silicate 
features . The details should be identical on a 
frequency scale (from Ref. 1) . 
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vs 

Moisture 


FIGURE 2. REFLECTANCE OF QUARTZ BEACH SAMPLES. The dry spectrum is 
qualitatively what one would expect (Ref. 9); the small 
dip at 8.6 pm is a combination band superimposed on the 
9.3 pm fundamental (taken from Ref. 8). 
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FIGURE 


O 



a) Spectra of tightly-bound 
water molecules. 


CHABAZITE 1938 COLO 



b) Spectra of loosely-bound 
or liquid water molecules 


COMPARISON OF SPECTRAL SHAPE OF WATER 
MOLECULES (FROM REF. 1). 
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= 1430 cm 




FIGURE 4. ACTUAL FORM OF NORMAL VIBRATIONS IN THE 
CO 3 —COMPLEX (FROM REF. 6 ). 


table ! 



General 


Observed 

Calculated 

Band 

position 

Assignment 

position 

position 

I 

2 55 m. 

I'l + 

2,5s/u-2.S0/( sh 

2.53 

II 

2 35jti 

3111 

235/i-2.30ms1i 

•2J3jtt 

111 

2 16 m 

Ul + 21*) + 

2,16;r-2.l2;<Sh 

2.MiU 



or 3i>, + 

2,16a('2.12/is>' 

2J3/( 

IV 

2 00/t 

2r, + ^ 

2.00M~UQfM 


V 

1 90 M 

I- 3^3 


I.S6/( 


(from Ref, 2) 
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APPENDIX B 

'calculation -OF SENSOR' NOISE LEVELS 

The physical consxderatxons leading to Eq. (4) and Eq. (5) in the 
main text are presented in this section. 

The reflected radxance xs given by; 

^ ^1 ^I ^ 

where is the bidijrectional reflectance (sr of the material 

-2 

and H is the sample irradiance (W cm ) . 

/ r 

The reflected radiant intensity is given by; 

Ir = A.cos(0j^) = A cos(Bp^) (W sr 

2 ' 
where A is the target area (cm,) and 

is the angle to the reflected ray from the norma] to 
the surface (radxans) 

Received power xs gxven by : 

A 

P(R) = (I^T^ + Ipath) (W) 

R 

2 

where A^ is the area of recexving optics (cm ) 

R is the range between receiver and target (cm) 
xs the atmospherxc transmission 
Ip^th path radxant intensity (W sr ^) and 

is the transmission of the optics 
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Rcos 



From the geometry illustrated above we see that 


Acos (0j^) A^ 


.2 


where f is the focal length of optical system (cm) and 

2 

A, is the area of the detector (cm ) 
d 


Hence, (reflectance change) - Ap^^^ 

R 


A ttD Tf 

o _ o _ 

2 2 2 

f^ 4f 4F^, 

N 


= Ap^ X Tf A '—r 
1 o A I 0^2 


AP^ = Ap^ 


T x,H_Aj 
IT o A I d 
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Since 


where 


we have 


Pi " — 
1 TT 


Pj IS the total diffuse reflectance 
d 


Ap_ • TT = AP , 
1 a 


T T H A 
H 


Hence, the ncTse equivalent power is given by: 


NEP = 


v-FaT 

a 


(i.e., Eq. 6 of Section 9) 


where B is the electrical bandwidth ^ 

and D is the detectivity of the detector- 

Since Ap- corresponds to a signal/noise ratio of one, the noise equiva- 
d 

lent reflectance, NEAp is given by 

/A, • B 4 

NEAP = 


^aVo^“' \lh 


(i.e., Eq. 5 of Section 7) 
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' APPENDIX C 

SIMULATION OF SENSOR NOISE THROUGH QUANTIZATION 
OF INBAND SIGNAL VALUES 

There are two, equivalent, ways in which the presence of a random 

-2 -1 

noise component of standard deviation of say 7 pwatts cm sr could 
have been simulated. 

The first way is to introduce a random noise signal of this magni- 
tude directly on the inband radiance values and. repeat the linear dis- 
criminant analysis for each atmosphere/ condition of observation several 
hundred times to determine an estimate of the average classification 
accuracy over the ensemble. We did not use this Monte Carlo technique 
because of the tremendous computational cost involved. 

The second way is based on the observation that what we are seeking 
is a method to establish a metric in the vector space of the inband radi- 
ance values such that one division along any axis corresponds to a noise 
level of a given probability. For example, at one-sigma noise level 
the units along each axis should correspond to the levels which can be 

discerned 68% of the time in the presence of Gaussian noise of say 
—2 —1 

7 ywatts cm sr 

■ ' -2 -1 
At first thought quantization into intervals of 7 ywatts cm sr 

appears to fulfill this role. However, since the inband radiance values 

are uniformly distributed over the quantization interval (i.e., any 

I ' 

value is equally likely) the standard deviation associated with this 
procedure is smaller than desired by a factor of /l2. 
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Specif xcally, the standard deviation, a, of a uniform distribution 
on the interval, (p — L/2, p + _L/2), is obtained from: 

U+L/2 

2 f 2 

a = 1/L • (x-u)-^ dx 

V-L/2 

y+L/2 
" y-L/2 




Hence, quantization of the inband radiance values into intervals of 

t 

width equal to t/12 • NEAL prior to performing the Linear discriminant 
analysis is an equivalent, but a much less costly way to predict sensor 
performance in. the presence of system noise. This procedure is the one 
we employed in this study. 
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